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Abstract 

We review some rigorous results on the equilibrium and non-equilibrium 
properties of superfluids and superconductors. 


1 Introduction and motivation 

The exact or rigorous study of equilibrium and non-equilibrium properties 
of superconductors has been initiated by G. L. Sewell in the 1980’s. We 
refer to chapter 7, part III, of his book [ Sew02 j for a beautiful introduction, 
and references to his own work on the subject, in particular jSew90j for the 
first proof of the Meissner effect in the bulk under the assumption of off- 
diagonal long-range order (ODLRO). In |SW09j the implications of ODLRO 
in superfluidity were explored, including spontaneous symmetry breaking 
(SSB). For a discussion of the relationship between ODLRO and SSB, see 
the very comprehensive book by Verbeure [Verllj . and for a recent study of 
the connection of those properties with Bogoliubov’s method of quasiaverages 
[bjTo] , see [WZ] , 

A recent, very stimulating historical review by Kadanoff |Kadl3| high¬ 
lights the coherence properties of the (’’slippery”) wave-functions in both the 
phenomena of superfluidity and superconductivity, a subject well-treated in 
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textbooks from the phenomenological point of view (see [ MR04] . chapter 5.2 
- the London theory). The microscopic basis of the London theory | LL35j is 
the property of ODLRO [Sew02j . We shall come back to this point briefly 
when dealing with superconductivity and the Meissner effect, but otherwise 
shall not discuss coherence properties any further. 

Superfluidity is the property of certain fluids, typically liquid Helium II, of 
flowing along pipes of extremely small diameter (ICR 2 cm or smaller) with no 
resistence (viscosity). Kadanoff, in a recent historical review |Kadl3] sharply 
questioned the relevance of Landau’s criterion to the superfluid property. The 
latter may be roughly stated in the following way: by the flow of a fluid along 
a pipe, momentum may be lost to the walls only if the modulus of the velocity 
|n| is greater than 

_ . e(p) 


where e(p) are the energies of the “elementary excitations” generated by 
friction. This indicates that the concept of superfluidity still lacks a clear 
and precise theoretical foundation. One reason, argued by Kadanoff, is that 
’’given the many mechanisms for broadening the distributions of both en¬ 
ergy and momentum, it seems very implausible that the Landau condition 
can begin to account for the very long-lived nature of the flow of superfluid 
Helium”. 

(Low temperature) superconductivity is the property of certain materi¬ 
als (e.g. mercury, tin, lead, aluminium) of becoming perfectly conducting. 
This is characterized by persistent currents, i.e., ideally everlasting cur¬ 
rents in superconducting rings, in the absence of external magnetic fields. 
Together with the phenomenon of persistent currents, the prototypical elec¬ 
trodynamic property of superconductors, now in the presence of an external 
magnetic field,, is the Meissner effect ( |MR04j . Chapter 5.2), one of the 
most spectacular effects in physics. It may be stated as follows: 

Meissner effect If a superconductor sample is submitted to a magnetic 
(applied) field H and then cooled to a temperature T below the transition 
temperature, there is a critical field H C {T ) > 0 such that, if \H\ < H C (T), the 
field is expelled from all points of the sample situated sufficiently far from 
from the surface, i.e., beyond a certain distance from the surface called the 
’’penetration depth”. 

Both (low temperature) superconductivity and superfluidity occur at very 
low temperatures, typically a few degrees Kelvin. 
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Much less understood than the coherence properties, and still a concep¬ 
tual challenge in several respects is the issue of whether the properties of 
superfluids and superconductors mentioned above, are equilibrium or non¬ 
equilibrium properties |Wrel5aj . This issue is also relevant to a third class of 
phenomena, that of magnetism: all three have in common the fact that 
they are macroscopic (i.e, due to infinite number of degrees of freedom) 
manifestations of quantum mechanics. Magnetism (diamagnetism) is also 
relevant to superconductivity, in the Meissner effect, and both ferro- and an¬ 
tiferromagnetism are relevant to certain aspects of superconductivity [Sew]. 
For non-equilibrium models of (disordered and ordered) ferromagnetism, see 
section 6.5.3 of [MW13j and references given there. 

Since the thesis of Niels Bohr (1912) and the work of van Leeuwen [Lee 2U> 
paramagnetism and diamagnetism were both known to require a quantum 
mechanical description as equilibrium phenomena, because the partition 
function has the property 


Z(A) = / dxdpexp(—/3H^(x,p) = Z(A = 6) 


if H A (x,p) = 


_ (p—eA ) 2 

2m 


+ V(x)i by a canonical transformation. 


We now ask: is the flow of a fluid in superfluidity or the persistent currents 
in superconductivity an equilibrium or a non-equilibrium phenomenon? 

In this review we should like to explain the exact results known about this 
problem. Since the natural framework for this study is the theory of infinite 
systems, which requires the use of algebras of operators, not well-known to a 
large portion of the community of theoreticians, we shall dwell mostly on the 
conceptually familiar finite systems, with the results on algebras of operators 
refered to and (hopefully) briefly clarified, without requiring from the reader 
any previous knowledge about them, except if they wish to master the details. 
The main focus will lie on the concepts, and their precise definition, which 
is a specially subtle matter in this connection. Certain equilibrium-non- 
equilibrium aspects (including the Sirugue-Winnink argument), and a simpler 
discussion of the application to the Girardeau model, are new. 

A crucial requirement assumed on the models treated here is Galilean 
covariance. For a different approach, see S14 . 
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2 General formalism, finite and infinite sys¬ 
tems 

We consider Bosons in translational motion (Fermions, although essential to 
superconductivity, will not need to be described explicitly, for reasons which 
will be explained later), to begin with in finite regions, which we take to 
be cubes, generically denoted by A. The thermodynamic limit will be taken 
along the sequence of cubes 

A n = [—nL, nL\ d of side L n = 2nL with n = 1, 2, 3, • • • (1-1) 

where d is the space dimension. The number L is arbitrary, and the sequence 
{n — 1, 2, 3, • • •} may be replaced by any subset of the set of positive integers. 
Let 

= -^per-(Aj) for J = 1,2,3,--- (1.2) 

denote the Hilbert space consisting of functions on R d , with / 6 L 2 (Aj) and 
such that f is periodic with period Lj = 2 jL with j = 1, 2, 3, • • • in each of 
the variables x = (x;) with i — 1, ■ ■ ■ , d, i.e., /(x* + 2 jL) = /(x*) with 1 = 
1, • • • , d and x = (x*), i = 1, • • • , d. Let 

J -a . denote the symmetrical Fock space over "H Aj with j = 1, 2, 3, • • • (1.3) 

We shall sometimes consider 

•Fkj,Nj = ® S) N'Hk j with j = 1, 2, 3, • • • (1.4) 

which denotes the symmetrized tensor product of 'H\ j corresponding to 
Nj particles in A j. Our local algebras of observables will be taken as the 
(von Neumann) algebras of bounded operators on Wy (or on Wy, jv, where 
Nj/\Aj\ = p, p is the density and |A| denotes the volume of A): 

Ay = £(Wy) (1-5) 

Such operators may be generated by the so-called Weyl operators, see jVerllj . 
By the above choice one has the isotony property 

C Aa; for k < / or A^ C A; (1.6) 

Our states will be assumed to be functionals of the local observables which 
on A a reduce to normal states, i.e., states of the form 

^a(A) = Tr u a (paA) VA e Aa , (1-7) 
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where p\ is a density matrix, i.e., a positive, normalized trace-class operator 
on "Ha; for the ground state, p A = |^a) (^a|, where |Ha) is a normalized 
vector in "Ha- A ground state will be denoted by 

UJ\ (Ha, ' Ha) ) 


and a temperature state 




Tr-H A (exp(-/3hT A )-) 

ZaW) 


/3>0, 


( 1 . 8 ) 


(1.9) 

where 

Z A (ff) = Tr HA exp (-/3H A ). 

Above, "Ha = Pa, si the (symmetrical) Fock space over L 2 (A). The generator 
of space translations on "Ha — the total momentum — will be denoted by 
Pa- Above, Ha is a ground state eigenvector of Pa on Ha: 

Ha^a = E a Q a , 


( 1 . 10 ) 

where 


Pa = inf spec(PA) • 


In addition, 


By thermodynamic stability, 


PaHa = 0. 


( 1 . 11 ) 

( 1 . 12 ) 


Pa > —c|A|, 


where |A| is the volume of A, and c is a positive constant. In general, E\ 
is of order of 0(— d|A|) for some d > 0, and in order to obtain a physical 
Hamiltonian satisfying positivity it is necessary to perform a renormalization 
(infinite in the thermodynamic limit) 

Ha H A = H\ — Pa . (1-13) 
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Define 


a\ )t (A) = exp(itH\)Aexp(—itH\) 

&k,x{A) = exp(Lr • PA)Aexp(—ix ■ P A ) 


(1.14) 

We now indicate how Galilean transformations act on finite systems. In A 
we consider a generic conservative system of N identical particles of mass m. 
In units in which h = m = 1, H A and P\ take the standard forms 

- y n a 

H a = ^ r + V(x l7 ...,x N ) (1.15) 

with V a suitable potential (satisfying certain conditions which will be spec¬ 
ified later) and 

N 

PA=-^Vr (1-16) 

r= 1 

with usual notations for the Laplacean A r and the gradient V r acting on the 
coordinates of the r-th particle. We assume that H\ and P\ are self-adjoint 
operators acting on = L 2 (A)f N , the symmetrized tensor product of 

"Ha corresponding to N particles, with domains D(H A ) and D(Pa), and 

D(P a )dD(H a ). (1.17) 

Let 

Si = {^f\ neZ d } 

and, given v G M d , let vh l ,l = ^h l ,l such that 

\kn L ,L ~v\= inf |k - v\ 
kesi 

and|fc Si)L | < |n| . 

If there is more than one kft Lt L satisfying the above, we pick any one of them. 
We shall refer to this briefly as the prescription. We have: 

lim Vft L) L = v, 

TV,L—4oo 


6 



where N, L —» oo will be always taken to mean the thermodynamic limit, 
whereby 


N —>■ oo, L —> oo, ~i = P with 0 < p < oo , 

where p is a hxed density. The unitary operator of Galilei transformations 
appropriate to velocity vh l ,l follows (upon restriction to the IV-particle sub¬ 
space of symmetric Fock space): 

U'l = exp (ivg L:L • (fi + ... + x N )) . 


(1.18) 

We shall assume that U'% maps D(H\) into D(H\). From now on we shall 
write v for vn L ,L- It follows that, on 

(uI)'h a u* = 

+ AEtf(A) , 




(1.19) 

where 

H\,v = + v ■ P\ 

(1.20) 

and 

N(v) 2 

A^(A) = 

(1.21) 


We shall need to describe infinite systems and their corresponding states. 
This necessity arises from the fact that 1) only in this limit a simple de¬ 
scription of the system is possible, independent of the peculiarities of finite 
systems (boundary conditions, etc.) 2) except in rare cases, for which surface 
terms become relevant even for very large systems (not met in the present 
article), this description is physically sound. 

For this purpose we define a so-called quasi-local agebra, which, in our 
case, due to the isotony property (1.6), is the limit as k —> oo, of *4.A fc in 
a suitable topology, which will be the weak topology, for reasons explained 
in (DS70] . We shall denote this algebra by A. The (inhnite-volume) states 
will be the thermodynamic limit of states (1.7), which will be denoted by oo, 
with a subscript /? or none, indicating whether we are dealing with thermal 
states (1.9) or ground states (1.8). They have the properties of positivity 
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(jj(A*A ) > 0 and normalization cu(l) = 1, inherited from (1.8) and (1.9), 
and are taken as defining properties of a state of an infinite system. An 
automorphism of an operator algebra is a one-to-one mapping of the algebra 
in itself, which preserves the algebraic structure. For the finite algebras the 
time- and space-translation automorphisms are given by (1.14). We assume 
that in the thermodynamic limit they give rise to automorphisms a t and <jg 
of A. 

We shall be interested in time-translation invariant states co(a t (A)) = 
ui{A) and homogeneous or space-translation invariant states: Vx G RAn o 
ag = co. 

A state co is associated to a triple (fi w , 7r^(kl), T-LJ) by the GNS construc¬ 
tion, see [ Sew86j . pg. 27, where tt u is a representation of A on a Hilbert 
space 

For time-translation and space-translation invariant states co , in terms of 
this representation, n bJ (at)(A)) = exp(itH ul )7T Ld (A) exp(— itH u ) for all A G A 
and similarly for the space translation automorphisms cxg(A) 7 T UJ (ag(A)) = 
exp(ix • P w )Aexp(—ix ■ P u ), where H^ and P u , the generators, are the (phys¬ 
ical) Hamiltonian and momentum, satisfying 


HMu, = o 

(1.22) 

PJX = o 

(1.23) 


Finally, we have that a t o a$ = otg o cc t , that is, the space and time 
translation automorphisms commute. We refer to |Sew86j . chapter 1, for a 
specially comprehensive introduction to the algebraic description of rnany- 
body (infinite) systems. 

After the renormalization (1.13), which is only possible if the finite- 
volume hamiltonian is bounded below (semibounded), we have H\ > 0, 
which leads in the thermodynamic limit to the condition 

H u > 0 (1.24) 

(1.22 also implies that co is a t - invariant. Concerning thermal states, the finite 
volume Gibbs states (1.9) satisfy certain analyticity properties known as the 
Kubo-Martin-Schwinger (KMS) boundary condition, see again (Sew86j for a 
comprehensive account. Although for a finite system and a given dynamics 
a\ it (see (1.14)), there is only one Gibbs state (1.9), this is not so for infinite 








systems, for which phase transitions may occur. It is therefore important to 
characterize equilibrium states of infinite systems by a more general condition 
[Sew86j . |Hug72| , [BR97] :let 0 < /3 < oo denote the inverse temperature. A 
state uj is a (cq,/3) - KMS state if, VA, £> e A, there exists a function Fa,b 
analytic inside the strip Dg = {A|0 < Tyz < /3}, bounded and continuous on 
its closure Vg, and satisfying the KMS boundary condition 

Fa,b( t) = u(Aa t (B)) and + i/3 ) = uu(a t (B)A) for all t G M. 


(1.25) 

A KMS state is cq-invariant. 

A state is a ground state if it is cq- invariant and satisfies the spectrum 
condition (1.24), and an equilibrium state at temperature (3 if it is a j3- KMS 
state. 

Although thermodynamic limits of Gibbs stes satisfy the KMS condition, 
it is not implied that such are the only ones. In fact, the KMS condition 
may itself be used to construct thermal states for infinite systems, as done 
by van Hemmen [ vH78j for BCS and related models of mean field type. 


3 A physical characterization of superfluid¬ 
ity and superconductivity 

The essential physical characterization of superfluid and superconductive 
states is the fact that they carry a current [Sew02] . One envisages, typi¬ 
cally, the following idealized situation. 

A superfluid or superconductor, initially in an equilibrium state (thermal 
or ground state), is set into motion by an external agent, e.g. a pressure 
gradient between two points of a pipe in the former case, a voltage gradient 
in the latter, which imparts to the system a steady macroscopic current, 
represented by a classical vector field: we refer to this situation as ’’the system 
in motion”. The main question we ask is whether the initial equilibrium state 
is also an equilibrium state of the system in motion. 

For simplicity, we consider first superfluids in a fixed inertial frame, the 
restframe of the system in the (equilibrium) ground state Ha- We assume 
that it is Galilean covariant, i.e., (1.19)-(1.21). The energy which must be 
imparted to the system to set it into motion with uniform velocity v is the 
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expectation value in Ha of 


UI) ] HkUI - h a 

— V • PA. + AE^ ;(A) 

which equals AE$(A). According to Landau, one should compare this en¬ 
ergy with the eigenvalues (or ’’elementary excitations”) of the Hamiltonian 
operator of the system in motion with velocity v, namely, U A y H A U A , which 
amounts to analyse the spectrum of H^ A , given by (1.20) : if it is positive 
for a given range of v, no dissipation occurs. For the relation of the above 
to the concepts of passivity and semipassivity, see the article by B. Kuckert 
[Kuc02j . 

By Galilean covariance, the momentum operator for the system in motion 
is given by 

Pj\.v = Pa + Afu (1.26) 

By (1.12), the system carries a current 

J A = Nv (1.27.1) 

which is the steady macroscopic current alluded to above. The system in 
motion is thus described by the Hamiltonian H A;i j given by (1.20). Due to 
(1.12), Ha is also an eigenstate of H A $ and therefore the former equilibrium 
state is also invariant under the dynamics of the system in motion. This fact 
generalizes to the thermal states (1.9), because H A $ commutes with H A , as 
well as to the thermodynamic limits of both ground and thermal states. 

The same considerations apply to the superconductor with a current, as 
remarked by Feynman in his lectures ( |Fey7 2|, 10.9). Instead of (1.27.1) we 
have 

J A = -Nev (1.27.2) 

where —e denotes the charge of an electron pair, and the excess energy, i.e, 
the energy furnished by the external agent (voltage) is AE$(A) ((10.56) and 
(10.57) of (Fby72| ). 

4 Nonequilibrium phenomena in superconduc¬ 
tivity and superfluidity 

We wish to inquire whether the initial equilibrium states of the system remain 
so in the presence of motion. For finite systems, this is not possible, both 
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for thermal states and for ground states. This does not, however, imply the 
same result for infinite systems. A Gibbs state at temperature /3 given by 
(1.9) cannot be a Gibbs state at the same temperature for the dynamics 
defined by (1.20), which we denote This follows from the fact that the 

requirement for all A £ Ak implies v = 0 by choosing 

A = exp(—/3H A ) — exp(— j3H A ^), which is possible by our choice of algebra, 
which is weakly closed. From the condition that limu;A(A) = limcnA^A), 
for all A £ A A) where lim denotes the thermodynamic limit, the conclusion 
does not follow, but we shall see shortly that, formulating the equilibrium 
condition as the KMS condition (1.25), the conclusion does follow. 

The reader may think that the above is a minor point of rigor, but perhaps 
the argument for ground states will convince him of the contrary. For finite 
systems, the ground state condition H\ > 0 does not imply the analogous 
condition H A ^ > 0 for the system in motion. Indeed, by (1.20), and (1.12), 

H A ,*Uf\n A ) = -l/2N\v\ 2 U^\Q A ) (1.27.3) 

(this argument is due to G.L.Sewell, see also [SW09] ). This is, however, 
not sufficient to show that the infinite system in motion has a physical 
Hamiltonian H w + v ■ which is not positive in the representation deter¬ 
mined by uj. (1.27.3) suggests that, ”in normal circumstances”, a has a 
spectrum contained in [— AE^(A), oo), which tends to the whole real line as 
N, L —» oo. In order to see what may go wrong, assume that the spectrum 
spec (Htf t a) = [— AE$(A), — Xn,l] U [0, oo), with X NjL —* oo as N, L —> oo. In 
this case, the negative part of the spectrum ’’disappears” in the thermody¬ 
namic limit, and the assertion that the spectrum of the physical Hamiltonian 
(for the infinite system) contains a non-empty set in the negative real axis 
does not follow. 

A further hint that (1.27.3) does not imply that the infinite system in 
motion has a physical Hamiltonian H u + v ■ P u which is not positive in the 
representation determined by to is given by the fact that, in the thermody¬ 
namic limit, the state G^IIIa) gives rise to a representation inequivalent (not 
unitarily equivalent, more precisely disjoint) to the one associated to uj. This 
is a well-known phenomenon, intrinsic to infinite number of degrees of free¬ 
dom, see chapter 1 of [ Sew86| or |MW13j . corollary 6.3, pg. 255: two different 
values of v, in particular v and 0, correspond to two different macroscopic 
densities, viz. momentum densities, by (1.26) 

It is therefore necessary to show that the spectrum of F/a ,v covers the 
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whole negative real axis in the thermodynamic limit, a fact which, as we 
shall see in the case of superfluidity, has a deep physical origin. 

In order to proceed, we consider both dynamics, one related to the inertial 
system (1.14) and the other to the system in motion. For a finite system, 
the latter’s corresponding automorphisms are given by 

a w (A) = e ‘tn,W Ae -UH,[K) ia^Ak, 


(1.28) 

and their infinite volume version will be denoted by cq^. For space transla¬ 
tions, by (1.26), the automorphisms ag are the same. 


4.1 A general approach to non-equilibrium states: the 
Sirugue-Winnink argument 

We assume that c o is time- and space-translation invariant, i.e., 


uj{a t (A)) = uj(A) 

(1.29) 

and 

u(a s (A)) = u(A) 

(1.30) 

where at and <7$ are commuting automorphisms of A Let 
of a Galilean covariant Boson or Fermion system, oj = (fl, ■ 

co denote a state 
D) and 

H U Q = 0 and = 0 

(1.31) 

Define the self-adjoint operators 


= H u T v • 

(1.32) 


for v G The corresponding automorphism group of A is the a ti g referred 
to above. We provide a sketch of proof of the following theorem, for those 
who only wish to understand the main argument. 

Theorem 1 Let u> = u>p satisfy the KMS condition with respect to a t for 
a given 0 < [3 < oo, where a t = a t g. Then, the quantum dynamical system 
(c 0 / 3 , a tt v, A ) defines a non equilibrium stationary state (NESS) whenever v ^ 
0 . 

Proof We must show that ujp does not satisfy the KMS condition with 
respect to the dynamics defined by a t $ . Sketch of proof: Assume the 
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contrary, i.e., that up satisfies the KMS condition with respect to cp^. Since 
it also satisfies the KMS condition with respect to cp by hypothesis, applying 
the KMS condition first w.r.t. cp and then w.r.t. cp^, we find, for A G A, 

Fa,B it) = up{a t {A)a t ${B)) = ujp(a tt ff(B)a t+i p(A)) 

i^t+i/3 i-^-)^t+i/3,viB ) ) 

= F AB (t + i/3 ). 


1.33) 

(1.33) leads straightforwardly to the fact that F a ,b is analytic in the whole 
complex plane and is therfore a constant. Varying A and B we gener¬ 
ate the whole Hilbert space and, from the definition of F, it follows that 
iUt^p)~ l U t ,u} P ^ = 1,, from which H u + v ■ P w = H u and, thus, v = 0. 

We now provide a complete proof, for those who wish to delve on the 
details. 

We actually obtain (1.33) only for A G A at , B G A ats , where A r is a 
norm-dense *-subalgebra of A consisting of entire analytic elements for r 
(see p3R87 j. Definition 2.5.20 and Proposition 2.5.22), The function Fa,b{z ) 
is, for A G A at , B G A ate , analytic in Vp = {z G C | 0 < < /3} and 

continuous on the closure T>p. By the three-line lemma (see, e.g., [BR97] . 
Proposition 5.3.5), it is uniformly bounded in Vp by ||A||||5||. Choosing 
sequences {A n } n > i, {B n } n > i with A n G A at , B n G A a± ^ such that 

\\An\\ < \\A\\- lim \\A n — H|| = 0 , 

n^oo 

|| B n || < || B ||; lim \\B n — B\\ = 0 , 

n—> oo 

one obtains that F A ,B n {z ) ->• F a ,b(z) uniformly in T>p (see [BR97J, Pg- 82), 
the limit function being therefore continuous and bounded on Vp, analytic 
in T>p, and satisfying the periodicity condition (1.33): 

Fa,b^) — Fa,b^ + i/3) Vt G M, V54, B G A . 

Furthermore, the following uniform bound holds as a consequence of the pre¬ 
vious inequalities , which themselves follow from Kaplansky’s density theo¬ 
rem, Theorem 2.4.16 of |BR87| : 

\F a , b (z)\ < ||A||||R|| MzeVp. (1.34) 
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Hence, by the Schwarz reflection principle, the function Fa,b extends uniquely 
to an analytic function in the whole of C, satisfying the bound (1.34), and is 
therefore a constant. From the definition (1.33), it follows that Fa.r{z) = 1 
for all ; £ C, whence the unitaries U t ^ e and U t)U p,v implementing the auto¬ 
morphisms a t and a t j in the GNS representation of ojp, satisfy 


( Ut,up ) Ut tU} p,v 


1 , 


upon using the cyclicity of f^. By our choices of a t ,a tj v : this implies that 
v = 0. q.e.d. 

Remark The above proof simplifies the original theorem in (SW70j . 

The above approach identifies what seems to be a central element in the 
issue equilibrium vs non-equilibrium in the present context: the same state 
and two different dynamics under which the state is assumed invariant (note 
that KMS states are necessarily invariant if the dynamics are implemented 
by an automorphism). The latter point accounts for the stationarity of the 
non-equilibrium states. 


4.2 The role of time reversal invariance: Bloch’s the¬ 
orem 

We now assume that the equilibrium state is time-reversal invariant, i.e., 
there exists an antiautomorphism 6 of the observable algebra such that 
u: o 9 = uj. Considering the representation tt u associated to uj by the GNS 
construction, assume that j(f) is a self-adjoint operator on a domain in 1-L W 
which includes fT, and which has the same property as the currents appear¬ 
ing in many-body quantum mechanics, i.e. Uej{f)Ug 1 = —j(f) where the 
antiunitary operatorf/g implements 6 in the GNS representation. Then we 
have 

Bloch’s theorem = 0 

Indeed, u(j(f)) = (u o 9)(j(f)) = u((9(j(f)))) = = 0. 

Theorem 1 generalizes Bloch’s theorem because it is an assertion on the 
state, not only restricted to the expectation value of the current, ft follows 
from both theorems that the persistent currents in superconductivity are 
definitely a non-equilibrium phenomenon. 


14 




4.3 Nonequilibrium nature of translational superflu¬ 
ids: the ground state 

Theorem 1 is valid for strictly positive temperatures, but not for the ground 
state. Moreover, the general basic assumption that the time and space trans¬ 
lations define automorphisms of the (weakly closed) quasilocal algebra, and 
even the existence of the physical Hamiltonian H w and physical momentum 
P u , very natural for any system, have been proved only for the free Bose 
gas and the BCS model |DS70j . It is therefore better |Wrel5a] to base one’s 
discussion on the Green’s function for the finite system [ BR97j . defined as 


G\^(A, B\ t, x) = uj A (Aa A ^ t (a AjS (B)) 

VA, B G Al and t G R 

It may be shown [BR97 ]. [Wrel5a] that there exist subsequences A nk such 
that the limits 


Gff(A, B\ t, x) = lim G An A A, B\ t, x) 

k K 

exist VA, B G A\/t G RVx G K d 

exist. Thus, a sequence A nfc such that the thermodynamic limit of the Green’s 
functions exists can always be found but is not unique: such nonuniqueness 
is related to the existence of phase transitions, which are expected for Boson 
systems at very low temperatures. 

We are now interested in considering the state a;, which is an equilibrium 
state at zero temperature (ground state), under the different dynamics, whose 
generator is H^ u . It is time-translation invariant, under this dynamics, but 
it is not necessarily an equilibrium state; precisely this situation has been 
studied before for KMS states. 

We have that |£)a) is an eigenstate of 77, 7 ,a of eigenvalue zero, so that 
stationarity is guaranteed and, indeed, follows naturally from the physical 
input without the need of performing a time-averaging on the state oj. 

We now define the energy-momentum spectrum ernsp A r k . i k of the finite 
system as the set of pairs (E rjA r kt L k (v),k Sj L k ) r ,s, where {E r ,N k ,L k (v)}r denotes 
the set of eigenvalues of 77y Ak = i7,7,v fc ,r fc , and {k s ^ k } s = S Ak the set of 
eigenvalues of the momentum PN k .L k ■ Above, when v occurs in the context 
of the finite system, the prescription is always understood. 
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Definition 1 emsp^v) is the set of limit points of emspN k) L k , possibly 
along subsequences. 

In |Wrel5a] we defined the following concept (see, e.g., [ BB03] , for the 
concept of support of a distribution): 


emspooiy) = The support of the Fourier transform of 
G$(A*, A; t, x) as a tempered distribution 

We also showed that it agrees with the support of the joint spectral family 
E(X,k) of in case these generators of time and space translations 

for the system in motion exist, and, furthermore, that 

emsp^iy) C emsp^iy) (1.35) 


This leads to 

Definition 2 emsp 00 (F) will be called the energy-momentum spectrum 
of the infinite system. If 3(A, k) e emsp^iy) such that A < 0, we shall say 
that the system describes a NESS. 

The last statement of definition 2 is justified by the fact that, when gener¬ 
ators exist, a zero-temperature state is a non-equlibrium state iff spec(H^ u )r\ 
(—oo,0) 7 ^ 0 by (1.24). Stationarity is included because, when generators 
exist, the corresponding state is time-translation invariant. 

Corollary 1 The corresponding ground state of the infinite system is a 
NESS If emsp^v) contains a point (A(/co),/co) with A(fco) < 0. 

4.4 Metastability: Landau superfluids as NESS 

Searching, now, for a metastability condition, i.e., a condition which might 
account for the long-livedness of the superfluid state, we are led to define 
subspaces of states ’’close” to the superfluid state, restricted to which H $ t a is 
expected to be positive. A natural proposal, connected to Kadanoff’s remark 
in the introduction related to the broadness in energy and momentum, is the 
following. Let 'H c I f L denote the subspace of T-Ln,l generated by of all linear 
combinations of vectors with simultaneous eigenvectors of 

Hn,l and Pn,l corresponding to eigenvalues E t) NjJv = 0 ) and kjy such that 

Ei,N,L < c for all i and \k h i\ < d 
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The restrictions in on the magnitude of the energy-momentum around 
the ground state above are due to the fact that we wish to express that 
there should be a sufficiently small energy-momentum transfer between the 
particles of the fluid and the surroundings (pipe). As remarked by Bayrn 
|Bay69| , the assumption that the energy of the excitations (in our case: the 
spectrum of Hfl) does not depend on the velocity of the walls is implicit. 
This assumption should be valid, but only for slow relative motion of the 
superfluid and walls. Accordingly, we pose: 

Definition 3 - A metastability condition The system defines a su¬ 
perfluid at T = 0 if 30 < v c < oo such that, whenever |u| < v c , there 
exists a subspace 1Z c fl l L of 'h% d L such that the eigenvalues {Ej NiL (v)'}j of the 

restriction of H NJj + vh l ,l • Pn,l = Hv.n,l to U% d L satisfy 
e c - d [y) = lim Ej NL (v ) > 0 for all j 

3 N,L-^oo ’ 

and that 

Wm, * {o} 

for some 0 < c < oo, 0 < d < oc. 

The last condition is a non-triviality requirement, which is expected to 
be satisfied in general (and is so for the Girardeau-Lieb-Liniger model). 

Definition 3 is related to Conjecture 2.2 of Cornean, Derezinski and Zin 
|DZ09] , ft may be expected to be a general characterization of metastability 
for superfluid Boson systems, independently of the presence or absence of 
BEC. An important issue is that definition 3 goes beyond thermodynamics, 
and for this reason Cornean et al define the inhmum of the energy spectrum 
(IES) in the thermodynamic limit e p (k), where p again denotes the density. 
We adopted a definition somewhat closer to the algebraic spirit, but both 
approaches have, of course, specific advantages. 


4.5 An illustrative example: the Girardeau-Lieb-Liniger 
model 


We start with the Girardeau-Lieb-Liniger model |Gir60j . |LL63] of N par¬ 
ticles (for simplicity odd) in one dimension with repulsive delta function 
interactions, whose formal Hamiltonian is given by (mass m = 1/2) 


N 

H n,l = ~Y^ 

i =1 


d 2 

(i dx ) 2 


N 

2 c( ) 5(xi — xfl with 0 < Xi, Xj < L . 

i,j=l 


(1.36) 
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The prime over the second sum indicates that the sum is confined to nearest 
neighbors. For the (standard) corresponding rigorous definition , see |Dor93j . 

The limit, as c —> oo, of equation (4), yields Girardeau’s model |Gir60j . 
It is the free particle Hamiltonian with Dirichlet b.c. on the lines x, = Xj, 
with quadratic form domain given in ( |Dor93j . pg. 353, (2.18)-(2.19)). It is 
straightforward that U'{ leaves this domain invariant and (1.17,1.19) hold. It 
is essential that we speak of the finite system, for which jv,l is bounded 
from below, and quadratic form methods apply; note that (1.17) is strict, 
because P\ is unbounded from below also for the finite system. 

For c —> oo, the b.c. on the wave-functions reduces to 

\l/(a;i,..., xn) — 0 if Xj — Xi, 1 <Xi,Xj<N, (1-37) 

and the (Bose) eigenfunctions satisfying equation (1.37) simplify to 

T b (xi, ..., x N ) = T f (xi, ..., x N )A(x 1 , ..., x N ) , (1.38) 

where dG is the Fermi wave-function for the free system of N particles con¬ 
fined to the region 0 < Xj < L, i — 1,..., N, with periodic b.c., and 

A{x i, ■ ■ ■, x n) = JJ sign(xj - x t ) . 

3<l 

Note that dG automatically satisfies equation (1.37) by the exclusion prin¬ 
ciple. Indicating the Bose and Fermi ground states by the subscript 0, it 
follows from equation (1.38) and the non-negativity of that 

*o =l*ol- 

Since A 2 = 1, the correspondence between T B and fiG given by equation 
(1.38) preserves all scalar products, and therefore the energy spectrum of 
the Bose system is the same as that of the free Fermi gas. Furthermore, 
dG is a Slater determinant of plane-wave functions labelled by wave vectors 
ki,i = 1 equally spaced over the range [— kp,kp], where kp is the 

Fermi momentum. Hence kp = which, in the thermodynamic limit 

reduces to 

kp = np. (1.39) 

The simplest excitation is obtained by moving a particle from kp to q > kp 
(or from —kp to q < —kp) thereby leaving a hole at kp (or —kp). This 
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excitation has momentum k = q — k F (or — (q — k F )), and energy e{k) = 
{k 2 - k 2 F )/2, i.e., 

e\k) = k 2 /2 + k F \k\. (1.40) 

This type of excitation must be supplemented by the umklapp excitations, 
which we consider in a more general form than [Lie63j . They consist in 
taking a particle from ( —k F — p ) to ( k F + q ) or ( k F — q ) to (—k F — p). We call 
the corresponding eigenvalues e 3 (k) and e 2 (fc) and consider just the latter in 
detail; for these 

2n 

0 < q < 2n(N — 1) and — < p 

and their momentum 

k = —2 k F — (p — q) (1.41) 

and their energies e 2 (/c) are given by 

2e 2 {k) = [(-k F - p) 2 - {k F - q ) 2 ] = [2 k F + (p - q)](p + q). (1.42) 

A different, equivalent choice involving holes may be made |Lie63] , 

We consider now a sufficiently large number of umklapp excitations; their 
momentum is given by (1.41) and thus opposite to v by our (arbitrary) 
choice (otherwise we would have taken the other set of umklapp excitations). 
Consider, now 

H V) n,l = H nl + vPn } l 

We now switch one particle from k F to — k F — p (corresponding to q — 0, 
p = ^4 in (1.41)). This corresponding elementary excitation aquires thereby 
the momentum 

, , 2vr 

k L = -2k F - — (1.43) 

jL/ 

Taking now a particle from k F — ^ to —k F — 2^4 (corresponding to q — 
^4, p = 2^ in (11.7), this again leads to the same elementary excitation 
momentum (1.43), and the resulting momentum of the two excitations is 
-2(2 k F ) + o(L). 

Performing now r successive umklapp excitations we obtain that the cor¬ 
responding excitation energy e(v, L ) of H Vi n,l = Hn,l + vPn,l is, by (1.43), 
with p = N/L, 

v r 

' Q 

e(v, L ) = (2 k F + 2tt/L) uJ^(-2 k F - 2n/L) = 

2—1 2 — 1 

= (2 k F + 27t / L) (27T / L -—^- 2tt / Lr) + v(—2 k F r — 2ir/L) 
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(1.44) 

which is —2kpvr + o(L), with kp given by (1.39). It is specially interesting 
that, for the finite system, (1.44) leads to the value 

Nv 2 

e(v,L) = -— + o(L) 

(notice that the mass m = 1/2), upon choosing r = [Lv/2tt\, that is, a 
sequence of umklapp excitations lead to the bottom of the spectrum (see 
(1.27.3)), up to small corrections. 

By further reasoning based on (1.40) |Wrel5a] , which we do not repeat 
here, we obtain: 

Proposition 1 For the Girardeau model, if 

0 < |u| < 2irp : 

the spectrum spec (H v> n,l) of H v N L contains the set £, defined by 
£ = {—2 kpvj + o(L)} with j = 1,2,3,--- 
Tlius, the quantities 

A j = —2kpvj with j — 1 , 2 ,... 

belong to emsp ^ in the sense of definition 1, which is ,thus, unbounded from 
below, and therefore the system describes a NESS by corollary 1. Further¬ 
more, adopting 

c = 2(7rp) 2 

and 

d = up , 

in definition 3, it follows that the system describes a superfluid at T = 0 in 
the sense described there. 

For c sufficiently large proposition 1 holds for the Lieb-Liniger model by 
continuity, but it is not straightforward to prove it for arbitrary c, although 
the result is expected, because a detailed analysis of the equations in the 
appendix of LL63j . not done there for the umklapp excitations, is required: 
there are, however, no problems of principle. 

Basically, with a linear behavior of the elementary excitations for \k\ 
small, c s \k\ +v-k > 0 for v sufficiently small. For large momenta and energies 
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other excitations spoil the picture; the Landau superfluid is a NESS. The 
umklapp excitations, which imparting large momentum and small energy 
to the particles, may be special to singular repulsive interactions in one¬ 
dimensional systems. Nevertheless, Seiringer and Yin SV()8i showed that the 
Lieb-Liniger model is a suitable limit of dilute Bosons in three dimensions, 
and, thus, a qualitatively similar picture may hold for more realistic systems. 
In fact, Lieb’s two branches of excitations have been seen experimentally, see 
[ TD03] . In general, for realistic systems, it is conjectured that the instability 
for translational superfluids is caused by rotons, which arise as local minima 
of the infimum of the excitation spectrum, see 1MR04| . pg. 246. These should 
also take place only for sufficiently large momentum, as a consequence of 
the properties of the liquid structure factor, if one assumes the Feynman 
variational wave function, see [ MR04] . pg. 256. For the Feynman variational 
wave function, see [DZ09j and |WdSJ05] 

The velocity of sound c s derived from the excitation spectrum c s = 
linifc^o in one dimension is related to the compressibility by c s = [(yy fy] - 
where P is the pressure P = -yyY For the Girardeau model (and also for 
the Lieb-Liniger model, but there the argument is less simple) infinite point 
repulsion makes E 0 proportional to which is proportional to 

N(N + l)(2iV + 1)/6L~ 2 . differentiating twice with respect to L we get 
-L'f L = o(£) which is nonzero in the thermodynamical limit. Comparing 
with free bosons, E$ = N/L 2 and the same calculation yields zero in the 
thermodynamic limit, or infinite compressibility. For this reason the free 
Bose gas is a bad model even in the limit of strong dilution. 

5 Equilibrium: the Meissner effect in super¬ 
conductivity 

As described in the introduction, the Meissner effect relates to supercon¬ 
ductors in the presence of an external magnetic field. The Hamiltonian is 
therefore not time-reversal invariant, and thus Bloch’s theorem is not applica¬ 
ble to the Meissner effect. In fact, strong arguments have been given (Sch60| 
that the latter is an equilibrium phenomenon. Of this remarkable review, 
|Sch60] , which appeared shortly before his untimely death in a plane crash, 
we also quote: ’’The most serious failing of the theory of BCS is, however, 
its failure to account for the electrodynamical properties of superconductors, 
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such as the Meissner effect and the persistent currents.” The reason for this 
is the BCS model’s lack of local gauge covariance, see |S ew02l , chapter 7, the 
introduction to |Wrel5b| and the footnote in |CS57j . 

The model we shall revisit was studied by Schafroth [ Sch55] as one of the 
very few locally gauge invariant systems possibly related to superconductiv¬ 
ity. Today, it may be viewed as a model for Schafroth pairs, which are known 
to occur in certain compounds , i.e., quasi-bound electron pairs localized in 
physical space, i.e., for which the spatial extension of the pair wave function, 
measured by the coherence length, is small compared with the average dis¬ 
tance between pairs. In this case all electrons of the band are paired, and 
the pairs form a dilute Bose gas. The present account from (Wrel5b| omits 
several technical details, e.g., regularity assumptions, and concentrates on 
just the essential points. The Hamiltonian may be written (e = 2eo, eo being 
the electron charge): 

H (A) = jjy + jA(x ) ) a -(x) ■ 

1(3 

•(V — yH(T))a(T)<£c 

(1.45) 

where a(x) and a*(x) are the basic destruction and creation operators on 
symmetric Fock space TT^Tf), with 77 = L 2 (K) the (one-particle) Hilbert 
space of square integrable wave functions on the cylinder K , assumed to be 
of radius R and height L, centered at the origin, of volume V = nR 2 L. 

77(A) is the (self-adjoint) second quantization of a one-particle operator 
on 77, with certain boundary conditions. We shall, however, use an extended 
version of this one-particle operator. The current density operator in this 
model is 

= j mom (x) + J Lon ^(x) (1.46.1) 

where 

jmom{x ) = : a*(x)Va(f) - a(£)Va(£) : (1.46.2) 

and 

Jlou,b ( f ) = - e —a*(x)a(x)A(x) (1.46.3) 

III 

Above, :: denotes the Wick or normal product. j Lon g{x) denotes the Lon¬ 
don part of the current. 
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Let (e p , eg, e z ) be a positively-oriented orthonormal basis adapted to cylin¬ 
drical coordinates. We assume that the magnetic induction B is uniform 
outside the cylinder, i.e., that the system is embedded in R 2 x [-L/2, L/2], 
and 

B(p, 6 , z) = Be z for p > R 

where B is a positive constant. The semiclassical (static) model will be 
defined by the above equation, together with the set of Maxwell equations 
for the magnetic induction B (using the MKS system): 

(V-R)(z) = 0 


(v x b){x) = no{n s ,j{x)n s ) 

where 

B(x) = (V x !)(£) 

and Qp is the ground state of H(A), assumed unique (see later). Note that 
the g.s. depnds itself on the solution to the Maxwell equation. 

At present no existence theorem for the semiclassical model is known, 
and we shall have to assume: 

Assumption 1 For any 0 < B < oo, the semiclassical model has a 
solution of the form 


f fee if p > R 
\ a(p)e e if p < R, 


The above (with K now replaced by the extended region) implies that 
the state |f2g) is of the form 






(A!) 1 / 2 


where |O 0 ) is the Fock vacuum (no-particle state), and 0 O denotes the nor¬ 
malized ground-state wave function of the one-particle operator 


H i = (P ~ eA(x)f 
a 2 m 

which, in correspondence to the embedding in R 2 x [—L/2, L/2] above, will 
be considered as an operator on 


U e = L\ R 2 , dp(p, 9 )) (8) L 2 ((—L/2, L/2)) 
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with 


dp(p, 6) = dp(p)dO with dp(p) = pdp 

Local gauge covariance is defined by 

U~ l H(A)U a = H(A + Vo) 

which is readily seen to be satisfied by the Hamiltonian H(A). The Maxwell 
equation is more precisely stated in the form (V x B)(f ) = ^ 0 (Og, j(/)Og); 
taking a sequence {/^} n > 1 , with /* —>■ <5(x) in the distributional sense, we 
may obtain the Maxwell equation in the form 

(V x B){x) = p 0 lim (%,/(/“)%) 

n—too 

in case the limit on the r.h.s. above exists, which is seen to hold. The 
standard choice of operators a, a*, implies a choice of phase a{x) = 0 or a 
real ground state wave function <f) 0 . Choosing a different a — a(x) instead, 
one obtains an additional term in the momentum density part of the current. 
For this reason, we should write, more precisely: 

e 2 

hon,^ = ~—a*(x)a(x)(A(x) - (Va)(x)) 

lit 

where Va is a gauge function which balances the gauge non-invariance of A, 
in order that a gauge invariant combination results. 

We assume that Neumann boundary conditions are imposed at z = ±L/2. 
h/T may be written 

(1-47-1) 

fcez 

where 

H ^ k) = ^-W-lk + 

\k + pa(p)\ 2 d 2 

p 2 dz 2 

(1.47.2) 




epB 


with a(p) given by 


if P > R , 






The above formula for the Hamiltonian leads us to expect that for B suf¬ 
ficiently small the ground state corresponds to k — 0, i.e., 0 O is independent 
of 6 , and that only the London part of the current contributes to the r.h.s. 
of the Maxwell equation. 

This is very difficult to prove, because a(p) also depends on this same 
wave function by the Maxwell equation. 

We thus pose: 

Assumption 2 There exists a constant Bo such that, for B < Bo, Hi, 
given by (1.47), has a ground state in the sector k — 0. 

Under this assumption the eigenvalue problem reduces to that of a radial 
Schroedinger equation with potential a(p) 2 and Maxwell’s equation becomes 
an integrodifferential for B(p), whose analysis by the contraction mapping 
principle in |Wrel5b| shows that it decays exponentially inside the cylin¬ 
der, with a ’’penetration depth” independent of the size of the sample and 
consistent with the physical data. 

This model is, therefore, not ’’exactly soluble”, as originally thought to 
be by Schafroth |Sch55j , on account of the fact that there is a ’’potential” 
cc(p) 2 , and, moreover, this potential is basically unknown! For this reason a 
proof of ODLRO for this model is very difficult. 

6 Conclusion 

We showed that current-carrying states cannot be thermal equilibrium states 
(theorem 1). In the case of superfluidity, it was demonstrated, together 
with results of |Wrel5aj . that ground states cannot be equilibrium states 
if a certain condition on the limit points (possibly along subsequences) of 
the energy-momentum spectrum of finite systems is met (corollary 1), and 
that in the latter case the systems describe non-equlibrium stationary states 
(NESS). In the Girardeau-Lieb-Liniger model it was shown that the condi¬ 
tions of corollary 1 are indeed met, and that the energy spectrum of the 
infinite system is unbounded from below (proposition 1). More important 
conceptually, this unboundedness from below was related to the existence of 
certain elementary excitations - the umklapp excitations - which derive from 
the existence of an effective ’’Fermi distribution” for the Boson system. The 
latter is, on the other hand, a direct consequence of the repulsive nature of 
the interactions, which, for realistic systems, is present in the van der Waals 
forces between neutral atoms. 
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The umklapp excitations, which imparting large momentum and small 
energy to the particles, may be special to singular repulsive interactions in 
one-dimensional systems. Nevertheless, Seiringer and Yin [ SY08j showed that 
the Lieb-Liniger model is a suitable limit of dilute Bosons in three dimensions, 
and, thus, a qualitatively similar picture may hold for more realistic systems. 
In fact, Lieb’s two branches of excitations have been seen experimentally, see 
[ TD03j . In general, for realistic systems, it is conjectured that the instability 
for translational superfluids is caused by rotons, which arise as local minima 
of the inhmum of the excitation spectrum, see [MR04 j. pg. 246. These should 
also take place only for sufficiently large momentum, as a consequence of 
the properties of the liquid structure factor, if one assumes the Feynman 
variational wave function, see | MR04j . pg. 256. For the Feynman variational 
wave function, see pDZ09] and |WdSJ05| . 

Another issue related to elementary excitations is that of metastability. 
The concept of local stability (of thermal states) introduced by Haag, Kastler 
and Trych-Pohlmeyer !IIKTP74 | is not suitable to describe stability of super¬ 
fluids and superconductors. This becomes specially clear in the deep anal¬ 
ysis of this concept in the case of ground states due to Bratelli, Kishimoto 
and Robinson [BKR78] . who showed that the condition of local stability in 
[HKTP74j when applied to the ground state is equivalent not to the ground 
state condition (1.24) (which would be the analogue of the equivalence con¬ 
dition proved - under certain assumptions - in [HKTP74 J in the thermal case, 
namely, the KMS condition (1.25)), but rather to the existence of a gap in 
the spectrum of the physical Hamiltonian H u . This is due to the possibility 
that the local perturbation creates an infinite number of infraparticles of 
infinitesimally small energy, a possibility which indeed occurs in superfluid¬ 
ity and superconductivity as a consequence of Goldstone Bosons arising from 
spontaneous symmetry breaking jVerllj . [SW09j . [WZj. At the same time, 
this suggests a clue to a possible definition of metastability, our definition 3, 
i.e., stability when restricted to a subspace spanned by ”a few” elementary 
excitations, essentially as suggested by Kadanoff [Kadl3] , see also [DZ09j . 
This is shown for the Girardeau-Lieb-Liniger model in proposition 1. In this 
connection, note that The ground state of a corresponds, by (1.27.3) , 
to the reversion of the velocities of a macroscopic number of particles. This 
explains the connection to metastability: very large momentum is necessary 
to connect the original ground state to this state. 

Lenard showed |Len64j that for the Girardeau model no BEC takes place. 
Therefore, superfluidity and BEC are conceptually different issues. However, 
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they may be related in the presence of Bose-Einstein condensation (BEC). 
The analysis of an interacting model with nonzero density in three dimensions 
exhibiting BEC and superfluidity remains as the most important, perennial 
open problem. 

For superconductors, there are no persistent currents in thermal equi¬ 
librium states by theorem 1, a generalization of Bloch’s theorem, which 
uses time-reversal invariance. Nevertheless, some of the best textbooks e.g. 
|Zim65] insist (pp.339-341) that persistent currents arise as a consequence of 
the gap in BCS theory (which is an equilibrium theory). The fact that Feyn¬ 
man’s Ansatz (1.27.2) requires Galilean covariance, which is broken in the 
BCS theory due to the truncation done to arrive at the model, also confirms 
that the BCS model is not adequate to analyse this issue. 

In spite of its lack of local gauge invariance, with serious consequences 
for the electrodynamics, analysed at length by Schafroth (Sch60j (see also 
jSew02] and |Wrel5bj ). the BCS model remains an important model in math¬ 
ematical physics, and the study of its mathematical structure, starting with 
the seminal work of R. Haag |Haa62] . has led to very important develop¬ 
ments, in the papers of Thirring and Wehrl [TW67 ] and Thirring [Thi68j . 
culminating in the comprehensive and more general approach of van Hemmen 
|ivH78j , who refers to Haag’s method as the ” Bogoliubov-Haag approach”. An 
important recent review [Id Si 61 summarizes extensive recent progress in the 
study of the mathematical properties of the BCS functional. 

Besides the open problem of persistent currents, which is a non-equilibrium 
effect, there are other effects which depend on the electrodynamics in super¬ 
conductivity and which are equilibrium effects, and, among these, the Meiss¬ 
ner effect plays a major role. Although Schafroth argues that the BCS Ansatz 
[CS57] in their approach to the Meissner effect is actually equivalent to the 
assumption of the existence of the effect itself, there is no general agreement 
on this issue, and, Kadanoff argued in his review |Kadl3] that the problem 
was solved by Anderson |And58] . who proposed that oscillations of the gap 
parameter produces extra states of the system (plasmons) which rescue the 
local gauge invariance of the theory. Although this may be true physically, 
it seems non-controversial that Anderson’s solution is not susceptible to rig¬ 
orous or exact bounds, and, in the sense of the present paper, the Meissner 
effect has really to be considered as an open problem, the more so that the 
BCS model in its original form is an exactly soluble mean-held model (see 
also j vH78 j and references given there). 

The model of Schafroth pairs introduced by Schafroth [Sch55 j to study 


27 






























the Meissner effect exactly, revisited in |Wrel5bj and reviewed in this paper 
seems to be the simplest, yet nontrivial enough to merit further study. We 
have already discussed why a proof of ODLRO for this models appears as 
(very) difficult. Assumptions 1 and 2 represent, therefore, the central open 
problem in this model: the existence of a unique solution of the London type, 
for sufficiently small values of the imposed magnetic held. 
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